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Abstract 

Hand-foot-and-mouth disease (HFMD) and herpangina are commonly prevalent illness in young children. They are 
similarly characterized by lesions on the skin and oral mucosa. Both diseases are associated with various enterovirus 
serotypes. In this study, enteroviruses from patients with these diseases in Korea in 2009 were isolated and 
analyzed. Demographic data for patients with HFMD and herpangina were compared and all enterovirus isolates 
were amplified in the VP1 region by reverse transcription-polymerase chain reaction and sequenced. Among the 
enterovirus isolates, prevalent agents were coxsackievirus A16 in HFMD and coxsackievirus A5 in herpangina. More 
prevalent months for HFMD were June (69.2%) and May (1 1.5%), and June (40.0%) and July (24.0%) for herpangina. 
Age prevalence of HFMD patients with enterovirus infection was 1 year (23.1%), 4 years (19.2%), and over 5 years 
(19.2%). However, the dominant age group of herpangina patients with enterovirus infection was 1 year (48.0%) 
followed by 2 years (28.0%). Comparison of pairwise VP1 nucleotide sequence alignment of all isolates within the 
same serotypes revealed high intra-type variation of CVA2 isolates (84.6-99.3% nucleotide identity). HFMD and 
herpangina showed differences in demographic data and serotypes of isolated enteroviruses, but there was no 
notable difference in amino acid sequences by clinical syndromes in multiple comparison of the partial VP1 gene 
sequence. 
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Introduction 

Hand-foot-and-mouth disease (HFMD) and herpangina, 
which commonly affect young children, are enterovirus 
infections causing a variety of exanthems. HFMD is a 
self-limiting exanthematous eruption characterized by 
vesicles in the oral cavity, mainly in the buccal mucosa 
and tongue, and peripherally distributed cutaneous 
lesions on the hands and feet. Herpangina produces 
multiple oral ulcers affecting predominantly the poster- 
ior part of the oral cavity only [1-3]. These diseases are 
associated with different strains of enteroviruses, such as 
coxsackievirus A (CVA) 2, 5, 6, 10, 16; coxsackievirus B 
(CVB) 1, 2, 5; and enterovirus (EV) 71 [4-8]. 
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Human enterovirus (HEV) genera containing the CVA, 
CVB, echovirus (ECV), and EV serotypes are transmitted 
mainly via the fecal-oral route and by contact with throat 
discharges or fluid from blisters [9]. Generally, HEV out- 
breaks peak during the summer and early fall, and various 
serotypes are often associated with a single outbreak [10]. 
Since 1993, when nationwide surveillance began in Korea, 
there have been reports of summer outbreaks of entero- 
viruses caused by ECV 5, 6, 7, 9, 13, 18, and 30; CVA 24; 
CVB 3 and 5; and EV 71 [11,12]. Especially, outbreaks of 
HFMD and herpangina caused by HEV infection were 
reported in 2009 in Korea [13,14]. 

Diagnosis of HEV infections is based on amplification of 
a highly conserved 5' non-coding region (NCR) that is 
widely-targeted in diagnostic procedures [15,16]. In 
addition to traditional virological methods to serotype 
HEV, reverse transcription-polymerase chain reaction 
(RT-PCR) based on amplification of the VP1 region have 
been recently developed [17-19]. Because the VP1 region 
is one of the main exposed regions of the viral capsid and 
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has been suggested to include a serotype specific antigenic 
neutralization site, the BC loop in this region has been 
implicated with viral antigenicity, and substitutions result- 
ing in conformational changes in this region are believed 
to play a role in host adaptation for HEVs [20,21]. 

In this study, to compare epidemic patterns of HFMD 
and herpangina, specimens from patients with HFMD 
and herpangina disease were collected and analyzed 
along with demographic data. Molecular detection by 5' 
NCR RT-PCR and sequencing of the VP1 region of 
HEVs were carried out. 

Materials and methods 

Diagnostic definitions 

Herpangina was defined as the presence of oral ulcers on 
the anterior tonsillar pillars, soft palate, buccal mucosa, or 
uvula. Patients with HFMD had oral ulcers on the tongue 
or the buccal mucosa, and vesicular rashes over the palms, 
soles, knees, or buttocks [22]. This study was conducted 
in accordance with ethical principles as formulated in the 
World Medical Association Declaration of Helsinki and 
approved by the institutional review board (IRB No. 
2012-48) of the Ethical Committee of Soonchunhyang 
University Cheonan Hospital. Additionally, informed con- 
sent was obtained from the parents on the patients' behalf 
who participated in the study and the parents of partici- 
pants also gave their consent to publish the data. 

Detection of enterovirus and molecular typing 

Stool specimens from 29 HFMD and 32 herpangina 
patients were collected from collaborative JeongGene 
Pediatrics hospital in Chungnam province, Korea during 
2009. The testing algorithm for detection and molecular 
typing of HEV has been previously described [12]. Viral 
RNA was extracted from the supernatant of infected 
cells using magnetic beads (Toyobo, Osaka, Japan). 
Extracted RNA was dissolved in 50 uL of nuclease-free 
water and stored at -70°C until used for RT-PCR. 

HEVs were assayed for in clinical samples using a 
AccuPower® Enterovirus Real-Time RT-PCR Kit (Bio- 
neer, Daejeon, Korea) based on 5' NCR of highly 



conserved region in the HEV genome, according to the 
manufacturers instructions [23]. If the sample was posi- 
tive, it was performed semi-nested PCR in VP1 coding 
region for molecular typing. Semi-nested PCR condition 
and primer sequences amplifying the VP1 coding region 
were described previously [18]. In the initial PCR, a 50 uL 
reaction mix containing 0.2 uM of primers 224 and 
222, 2 U of Taq DNA polymerase (Promega, Madison, 
WI), 100 uM concentrations of mixture of dNTPs, 
and 2 uM MgCl2 was amplified by 40 cycles of 95°C for 
30 sec, 42°C for 30 sec, and 60°C for 45 s. One microliter 
of the first PCR product was added to a second PCR for 
semi-nested amplification. Fifty microliters of a reaction 
mix containing 0.2 uM of primers AN89 and AN88, 2.5 U 
of Taq DNA polymerase (Promega), 100 uM concen- 
trations of a mixture of dNTPs, and 2 uM MgCl2 
was incubated at 95°C for 6 min prior to 40 amplifi- 
cation cycles of 95°C for 30 s, 60°C for 20 s, and 72°C for 
15 s. PCR products were purified using the QIA quick 
PCR purification kit (Qiagen, Valencia, CA). Purified DNA 
was added in a reaction mixture containing 2 uL of Big 
Dye terminator reaction mix (Applied Biosystems, Foster 
City, CA) and 2 pmoles of AN88 and AN89 primers. Se- 
quencing reactions were subjected to initial denaturation 
at 96°C for 1 min and 25 cycles consisting of 96°C for 
10 s, 50°C for 5 s, and 60°C for 4 min in a Gene Amp PCR 
system 2700 (Applied Biosystems, Foster City, CA). The 
products were purified by precipitation with 100% cold 
ethanol and 3 M sodium- acetate (pH 5.8), and then loaded 
on a model 3100 automated genetic analyzer (Applied 
Biosystems). The molecular type of each isolates was 
determined by the serotype of the highest scoring strain in 
Genbank using the Basic Local Alignment Search Tool 
(BLAST); that is, the sequence of the HEV strain that gave 
the highest nucleotide similarity value with the query se- 
quence [24]. 



Sequence analysis 

Sequence analysis was carried out for five strains of 
CVA 2, 16 strains of CVA 5, 18 strains of CVA 16, two 



Table 1 Numbers of enterovirus serotype isolated from patients with HFMD and herpangina 



Serotypes 



HFMD 



Herpangina 





Isolate numbers 


Percentage of subtotal 


Isolate numbers 


Percentage of subtotal 


Coxsakievirus A2 


0 




5 


20.0 


Coxsakievirus A5 


3 


11.5 


13 


52.0 


Coxsakievirus A16 


16 


61.5 


2 


8.0 


Coxsakievirus B1 


1 


3.8 


1 


4.0 


Enterovirus 71 


4 


15.4 


2 


8.0 


Untypable 


2 


7.7 


2 


8.0 


Total 


26 


100 


25 


100 
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strains of CVB 1, and six strains of EV71 (Table 1). Nu- 
cleotide and deduced amino acid sequences of all HEV 
isolates were compared with the reference strains using 
CLUSTAL W (version 1.81) and Megalign (DNASTAR) 
[25]. The phylogenetic relationships among the VP1 
sequences of each virus isolate were inferred using 
MEGA software v. 5.05. Maximum Composite Likeli- 
hood was used as the substitution method, while the 
neighbor- joining method was used to reconstruct the 
phylogenetic tree [26]. The reliability of the phylogenetic 
tree was determined by bootstrap re-sampling of 1,000 
replicates. 

Nucleotide sequence accession numbers 

The HEV sequences reported here were deposited in the 
GenBank sequence database under accession numbers 
JF773151 to JF773197. 

Results 

Enterovirus detection and molecular typing 

Among stool specimens from 29 HFMD and 32 herpan- 
gina patients, the approach detected 26 HEV isolates 
(89.7%) in HFMD and 25 HEV isolates (78.1%) in her- 
pangina. The types in the 26 isolates from HFMD 
patients were mainly identified as CVA 16 (16 isolates, 
61.5%), EV 71 (four isolates, 15.4%), CVA 5 (three iso- 
lates, 11.5%). The types in the 25 isolates from herpan- 
gina patients were mainly identified as CVA 5 (13 
isolates, 52.0%), CVA 2 (five isolates, 20.0%) (Table 1). 

Epidemiological features 

In temporal distribution, HEV positive samples from 
HFMD patients comprised three isolates (11.5%) in May, 
18 isolates (69.2%) in June, one isolate (3.8%) in July, 
two isolates (7.7%) in August, one isolate (3.8%) in 
September, and one isolate (3.8%) in October. The HEV 
positive samples from herpangina patients comprised 
five isolates (20.0%) in May, 10 isolates (40.0%) in June, 
six isolates (24.0%) in July, three isolates (12.0%) in 
August, and one isolate (4.0%) in September (Figure 1). 
Concerning the age of HFMD patients with enterovirus 
infection, two isolates (7.7%) were recovered from 
patients <l-year-of-age, six isolates (23.1%) at 1-year-of- 
age, four isolates (15.4%) at 2-years-of-age, four isolates 
(15.4%) at 3-years-of-age, five isolates (19.2%) at 4-years- 
of-age, and five isolates (19.2%) at >5-years-of-age. 
Concerning the age of herpangina patients with en- 
terovirus infection, two isolates (8.0%) were recovered 
at <l-year-of-age, 12 isolates (48.0%) at 1 -year- of- age, 
seven isolates (28.0%) at 2-years-of-age, three isolates 
(12.0%) at 3-years-of-age, and one isolate (4.0%) at 
4-years-of-age (Figure 2). 



Number 


20 






— ♦— HFMD 








A -■— Herpangina 


16 


/ \ 


14 


/ \ 


12 




10 




8 




6 




4 




2 




0 






May Jun Jul Aug Sep Oct 


Figure 1 Temporal distribution of enteroviruses isolated from 


patients with HFMD and herpangina. 



Sequence analysis of enteroviruses 

Based on the VP1 amino acid sequence comparison for 
all isolates, all sequences represented the same intra- 
serotype, except for the Kor09-CVA5-309cn (herpan- 
gina) isolate, regardless of syndrome (Figure 3). The 
pairwise VP1 nucleotide sequence alignment of all iso- 
lates was compared within the same serotypes; CVA2 
isolates displayed 84.6-99.3% nucleotide identity within 
the intra-type. CVA5, 16, CVB1, and EV71 isolates 
showed 98.7-100%, 89.6-100%, 99.4%, and 98.9-100% 
identity, respectively. The phylogenetic relationships of 
Korean isolates and reference strains are shown in 
Figure 4. 

Discussion 

HFMD and herpangina were epidemic in Korea in 2009 
[13,14] and are commonly prevalent illness in young 
children. They are similarly characterized by lesions on 
the skin and oral mucosa. Also, these two diseases are 
associated with various enterovirus serotypes. In HFMD, 
it has been reported that the major causative HEV types 
are CVA 16 and EV71; these types were also more preva- 
lent serotypes in Korea in 2009 [27,28]. Generally, CVA6 
and CVA 10 are the more prevalent serotypes in herpan- 
gina; however, CVA5 and CVA2 were the prevalent her- 
pangina serotypes in Korea in 2009 [29]. In multiple 
comparison of partial VP1 gene in the intratypes, no 
notable variation in amino acid sequences by clinical 
syndrome was evident. The VP1 region containing the 
BC loop is associated with viral antigenicity and substi- 
tutions of amino acid sequences in this region are 
believed to be important in host infection of HEVs 
[20,21,30]. No significant clustered relationship was 
found according to the clinical syndromes. When the 
pairwise VP1 nucleotide sequence alignment of all iso- 
lates was compared within the same serotypes, CVA2 
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Figure 2 Age distribution of enteroviruses isolated from patients with HFMD and herpangina. 



isolates showed an appreciable variation in nucleotide 
identity (84.6-99.3%) within the intra-type. Deduced 
amino acid sequences of all strains were in the same 
intra-serotype, except for one isolate, regardless of 
syndromes. 

There was no significant difference in age distribution be- 
tween herpengina (3.74 ± 2.48 years) and HFMD (3.54 ± 
2.05 years) in a previous report [31]. But, in this study, age 
distribution of HEV positive in herpangina patients was 
relatively lower than in HFMD patients. In the 2009 Korean 



epidemic, enterovirus infections were most prevalent in July 
[32]. However, temporal distribution of HEV positive in 
herpangina and HFMD peaked in June, with herpangina 
gradually decreased thereafter while HFMD decreased 
rapidly. 

In conclusion, HFMD and herpangina showed diffe- 
rences in demographic data and serotypes of isolated 
HEVs, but there was no notable difference in amino acid 
sequences by clinical syndromes in multiple comparison 
of the partial VP1 gene sequence. 
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Figure 3 Comparison of the deduced amino acid sequences of VP1 region of 47 enteroviruses isolated from patients with HFMD and 
herpangina. The rectangular parts are BC loop previously reported. 
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Figure 4 Phylogenetic analysis based on VP1 region of enteroviruses isolated from patients with HFMD and herpangina. 

Nucleotide sequences were analyzed by the neighbor-joining method. The numbers at the branches indicate the bootstrap values for 1,000 
replicates. 
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